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У вивченні геодиначічних явищ земної кори за даними гравіметрії, як 
одного з найбільш доступних засобів, доводиться зважати на дві обставини: 
завдяки | технологічному прогресу (удосконаленню апаратури 1 
топогеодезичного забезпечення СР5) зросла точність 1 продуктивність 
гравіметричних знімань, а через економічний спад різко знизився обсяг 
вимірювань. 

Перша обставина сприяє перегляду методів опрацювання отриманих 
даних, зокрема, точнішому обліку поправок Буге | 1, а друга, за необхідності 
виявляти глибинні джерела аномалій, потребує перегляду методики 
вимірювань, щоб урахувати тонкі особливості (варіації) гравіаномалій. Ці 
особливості можна задіяти через додаткову змінну - час. 

Світові "тренди" геофізичних спостережень тяжіють до 4) моніторингу 
|2| досліджуваної площі: вивчають еволюцію гравітаційного поля за час 
експлуатації ділянки чи в інтервалі її динамічної активізації; деякі аспекти 
моніторингу наведено в |3|. Застосування повторних спостережень у певній 
галузі для гравіметрії не нове, хоча безперервними в часі їх назвати важко, на 
відміну від сейсмометрії. 

У разі проведення високоточних вимірювань, параметром яких є сила 
тяжіння, варто обов'язково досліджувати й ураховувати її варіації (про що 
зазначено ще в |4|), щоб увести відповідні поправки у гравітаційні знімання 
різних років. Зокрема, виявлені неприпливні квазіперіодичні варіації (КІІВ) 
сили тяжіння на новій методологічній, методичній 1 метрологічній основі Й 
обгрунтовано їхню природу |4|. Внаслідок кореляції великомасштабних 
густинних неоднорідностей мантії 1 припливних параметрів Землі |5| 
можлива крос-кореляція КІІВ 1 цих параметрів, особливості якої ще потрібно 


вивчати. 
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Крім того, виявлено |6| неприпливні варіації вертикального градієнта 
(різниці відліків значень сили тяжіння на різних рівнях зі сталим перепадом 
висот, далі - ВВГ) сили тяжіння. На перший погляд, вони мають логічно 
випливати з варіацій сили тяжіння, проте розрахункова амплітуда таких 
варіацій дуже мала, тоді як реальні амплітуди на Кілька порядків вищі. 
Основні періоди неприпливних варіацій сили тяжіння, як 1 ВВГ, корелюють 
з періодами власних коливань Землі. Довготривала складова цих варіацій 
(зарік) корелює з положенням осі обертання Землі відносно Сонця 
(мінімальні амплітуди - взимку, максимальні - у квітні-травні 1 ЖОВТНІ- 
листопаді). Хаотичні флуктуації руху Сонця 1 Місяця можуть породжувати 
неприпливні ВВГ, однак доведення цього припущення - справа майбутнього. 

Зв'язок аномальних змін ВВГ із проявами землетрусів у вигляді 
мінімуму ВВГ задовго до початку поштовхів дає підстави припустити, що 
джерелом власних коливань Землі є не землетруси, а навпаки, власні 
коливання є спусковим гачком землетрусів. Це дає змогу простежити 
динаміку земної кори за критерієм середньоквадратичної похибки зміни 
амплітуди одиничного вимірювань ВВГ. Можливість тлумачення аномалій 
ВВГ як провісників землетрусів ще потребує опрацювання. 

Виділяють корисний сигнал шляхом усереднення у Кковзному вікні, 
спектральним аналізом (періодограми", оцінка спектральної густини). Крива 
зміни ВВГ - це складні квазіперіодичні флуктуації з різними періодами, 
амплітудами 1 фазами (рис. 1), 1 виділити певну залежність за варіаціями з 
окремих кривих неможливо. Однак, підсумовуючи криві варіацій за тривалий 
відтинок часу (від тижня до року), можна виділити певні періодичні варіації: 
з кривих ВВГ У,, знімають для кожного хвилинного відтинку часу значення 
ординат 1 для кожної хвилини 1-годинного інтервалу (збігається з основним 
фоном власних коливань Землі 54-60 хв) їх усереднюють, що взаємо- 
компенсує випадкові варіації чи ті, що мають інший період. 

У праці |6| запропоновано задіяти середньоперіодні ВВГ для визначення 
добротності мантії Землі: зі зменшенням періоду коливань на вимірювання 


впливають ближчі до поверхні шари геологічного розрізу, що позначається 


Ї Коефіцієнти у розкладенні часового ряду за гармоніками ехр Фгісі/ по , що виражає 


амплітуду певної гармоніки. 
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на картині ВВГ. У верхній частині розрізу ця картина ускладнена через 
більшу мінливість пружних властивостей гірських порід, зумовлену зміною 
літології 1 цілісності шарів. Короткоперіодні ВВГ неможливо виміряти через 
те, що крок дискретизації спостережень часто більший від періоду коливань: 
за такого співвідношення ці коливання зливаються з випадковими завадами 


вимірювань. 
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Рис. 1. Графіки варіацій У,,;: прямі та усереднені за місяць; ромбами позначено 


катастрофічні землетруси | 61. 
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Зв'язок дисперсії ВВГ з розломними зонами можна використати для 
оцінки тектонічної активності розломних зон (регіонального 1 локального 
(нженерно-будівельного) сейсмотектонічного районування): інтенсивність 
похибки залежить від активності розломів, якщо ці зони "заліковані" (масив 
стає майже однорідним), то така зона не виявиться у дисперсії ВВГ. 
З урахуванням того, що сила тяжіння є функцією маси, а вертикальний 
градієнт - функцією зміни густини, отримаємо нові можливості для вивчення 


глибинної будови Землі. 
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Рис. 2. Одинадцяти-денний ряд вимірювань сили тяжіння: без фільтрації, зі 


смуговою фільтрацією 1 з поправкою на припливний ефект. 


Урахування значень КІВ та ВВГ сили тяжіння є неодмінною передумо- 
вою для застосування гравітаційного моніторингу, про що не зазначено у | 21. 
Вивчення поведінки цих величин доцільно розпочати на профілях, наведених 
у 141. 

Гравітаційний моніторинг - це ряд періодично повторюваних у реаль- 
ному часі безперервних протягом фіксованого часового відтинку (рис. 2) 


мікрогравіметричних вимірювань гравітаційного поля й опрацювання його 
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даних з урахуванням впливу навколишнього середовища та сфери 
застосувань. Значення часового відтинку залежить від якості вимірювань, 
міри невизначеності результатів спостережень, динаміки (амплітуди 1 
частоти) поля. 

Його фізична основа - безперервний зв'язок динаміки гравітаційного 
поля 1 параметрів середовища. Просторовий розподіл варіацій вертикальної 
похідної У, потенціалу сили тяжіння прямо корелює з площинним 


розподілом густин, а часові варіації У, чітко визначають вертикальні варіації 


насичення флюїдів. 

Апаратна основа моніторингу - спільні великомасштабні вимірювання 
перевищень рельєфу за даними СР5 1 абсолютних значень сили ТЯжінНнняЯ. 
Відносні вимірювання сили тяжіння дешевші, проте обмежені прив'язкою до 
опорної мережі 1 синхронністю обліку "сповзання нуля". У наших умовах, 
однак, такі вимірювання за належного методичного забезпечення |4| мають 
перспективу. 

Мета моніторингу -- оцінка глибини до джерела аномалій 1 зміни об'єму за 
даними деформації рельєфу (з даних  СР5) - потребує знання 
приповерхневого розподілу мас (з даних гравіметрії). 

У дослідженнях з 4Р гравімоніторингу" визначена амплітуда сигналу в 
діапазоні від 20 до 30 мГал, у межах діючих вулканів вона зростає до 300- 
600 мГал, а в межах водних басейнів - 200-250 мГал у разі нелінійного 
обліку "сповзання нуля". 

Виділяють слабкий сигнал на тлі перешкод за відомими у гравіметрії 
методами кореляційного аналізу й обчислення складових градієнтів у прямих 
вимірюваннях сили тяжіння. У |2| виділяють корисний сигнал з різницевої 


аномалії Ля З Дб 7 УП - Д8дер - Д8у» де ДА - різниця між сусідніми 


відліками в часі; у - 308,6 мГал/м - поправка у вільному повітрі; Й, - 


вертикальний зсув (відносне підняття опускання); ДЄде; - аномалія Буге 


деформації (зміни об'єму через стиск середовища навколо джерела, що 


означає зсув густинних меж У неоднорідному середовищі); 


2 Усі статті спецвипуску СеорПузіс8. 2008. Моі73. М. 6 присвячені гравітаційному 
моніторингу. 
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ДЯ у - 2ЯСРуФО, «42ф0. - вплив грунтових вод. Також порівнюють 


результати фільтрації з даними, отриманими з контрольного пункту, близько 
розташованого від ділянки спостережень. 

На рівнинних ділянках на території України корисний сигнал матиме 
малий градієнт, тому треба враховувати вплив перешкод іншими методами 
(зміна геометрії мережі спостережень, часова фільтрація, обчислення 
похідних тощо), запропонованими у |4|. Вивчення метрологічних нюансів 
доцільне на базі геодинамічного полігона з готовою інфраструктурою 
(можливо, на базі однієї з ділянок Карпатського відділення ПФ). 

Перш ніж перейти до умов стійкого тлумачення даних, зазначимо, що 
вимірювання на регулярній мережі й перерахунок значень за інтегралом 
Пуассона придатні для регіональних досліджень |З|, а в локальних умовах 
гравімоніторингу |9| мають низку недоліків | 10|. Іноді організація регулярної 
мережі неможлива, а перерахунок з нерегулярної мережі на регулярну 
складніший, ніж обернене відновлення будови середовища. 

Розв'язування обернених задач гравіметрії з даними на псевдорегулярних 
мережах із застосуванням моделей середовища типу "безкінечний профіль" 
утворює погано зумовлені системи лінійних рівнянь із беззмістовними 
результатами. З огляду на це 1 на іншу Й коротку довжину реальних профілів 
доцільний альтернативний підхід |10| - система лінійних інтегральних 


рівнянь зі швидкоспадними ядрами: 


-Ї 
зем в і Т9ОХ дез є 
а «г 2 п Є. 
сов 250 Є 350 Є губ 2 соб 250 Є 350 во по оо 

З урахуванням цього методика |9| дієва за умов, що регіональний фон -- 
поліном першого степеня, відомі густини 1 положення меж ТІіл, що тяжіють, 
на поверхні, а тіла близькі або мають спільні контакти. Однак ми пропонуємо 
замінити згадані у методиці |9| програми на комплекс | 1 1 | 1 програми з (10). 

У |2| не враховано й таку особливість варіацій гравітаційного поля: 
величини його похідних залежні від флуктуацій аномальних атмосферних 


мас, снігових мас, рівня грунтових вод, лісистості 1 змін рельєфу через 
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техногенну діяльність | 12|. Урахувати ці ефекти можна, вносячи відповідні 
поправки |12| у розв'язок прямих задач гравіметрії на дослідних територіях. 

Необгрунтоване спрощення моделей середовища для зниження 
неоднозначності у багатьох випадках є причиною неправильних обчислень 
геометрії джерел, вертикального 1 латерального розподілу густинних 
неоднорідностей,  собливо у | випадках, коли середовище навколо 
аномального джерела далеке від припущень про однорідність. Надійне 
тлумачення одержують за відомої геометрії тіл, що тяжіють (ІЗ 
сейсморозвідки), 1 комплексної інтерпретації поля сили тяжіння та 
деформацій рельєфу. 

Враховувати наявну апріорну інформацію про середовище пропонуємо 
шляхом: 

Г) побудови належних модельних зображень (зіркові області відомої 
густини усередині компактних множин у банаховому просторі даних) | 101; 

2) додання у функціонали типу відхилу у регуляризувальних алгоритмах 
стабілізаторів диференціального вигляду, власні функції яких збігаються з 
власними функціями вихідних операторів (131. 

Розв'язки обернених задач за даними гравітаційного моніторингу 
доцільно отримувати за регуляризувальними алгоритмами з урахуванням 
розробок |10, 11|. З них можна відновити загальну картину варіацій густини, 
проте не абсолютні значення густин. Окрім завдань геодинаміки, ці розв'язки 


можна залучити для моніторингу стану екосистем. 
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Рие го Фе ресиПагіпез ої єгамту пттопіїогіпе 165 Чака ргосе85іпе оп Ше 5Погі ргоНіе5 18 
оНегей го ехесиїе Бу Ше ехрге8510п5 уліб Га5і десаутє, соге5, апа єгамту Писішацопя 
де го трасі ої 5пааї іпіепзтуе єеоррувбісаї! Гасіог5 опе 5пошіа їаКе їпіо ассоцпі Бу Фе 


РмиПі8 іесрпідцез. 
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ОСОБЛИВОСТІ ПРОСТОРОВО-ЧАСОВОЇ СТРУКТУРИ 
ВІКОВИХ ВАРІАЦІЙ МАГНІТНОГО ПОЛЯ ЗЕМЛІ 
НАПРИКІНЦІ ХХ - НА ПОЧАТКУ ХХІ СТ. НА ТЕРИТОРІЇ ЄВРОПИ 


В. Максимчук, Ю. Городиський, Д. Марченко 


Карпатське відділення Інституту геофізики ім. С.І.Субботіна НАН України 
Е-пааї! упаакзутісрик Фсб-ієрб.Їміу. па, уогеогод Фсбі-єрі.Туту ца 


Особливістю магнітного поля Землі є його змінність у часі. Варіації 
магнітного поля мають як внутрішні, так зовнішні (тіоносферні) джерела. 
Розглянемо вікові зміни, або віковий хід (5У), геомагнітного поля 
внутрішнього походження. Дані 0 про | просторово-часову структуру 
геомагнітного поля необхідно враховувати під час проведення магнітних 
знімань різного призначення, побудови карт аномальної частини магнітного 
поля, у навігації тощо. Головним джерелом інформації про просторово- 
часову структуру вікового ходу є магнітні обсерваторії (МО) 1 дані 
вимірювань на пунктах вікового ходу. Побудовані за цими даними карти 
вікового ходу компонент геомагнітного поля (карти ізопор) дають уявлення 
про морфологію 5У в досліджуваних регіонах, а їхнє порівняння для різних 
епох дає змогу виявити центри (фокуси) вікового ходу та вивчати їхню 
кінематику. Найліпші умови для вивчення 5У є на території Європи завдяки 
доволі щільній мережі МО. 

Просторову структуру вікового ходу на території Європи досліджено в 
П1-3|. У цих працях з'ясовано, що у другій половині ХХ ст. в Європі 
простежувалась доволі динамічна структура ЗУ: у 60-х роках відбувся розпад 
Каспійського (ТГранського) фокуса 5У 1 формування в 70-х роках нового 
фокуса 5У в Центральній Європі | 1 |. 

Наша мета - побудова карт вікового ходу геомагнітного поля для 
території Європи наприкінці ХХ та на початку ХХІ ст. 1 вивчення на ЦІЙ 
підставі особливостей його просторово-часової структури. 

Даними для побудови карт ЗУ є середньорічні значення компонент 
геомагнітного поля на МО Європи за період 1985-2005 рр. |4|. Карти 5М 
були побудовані для всіх компонент поля чотирьох 5-річних інтервалів з 
1985 по 2005 рр. (З огляду на обмеження обсягу повідомлення ми наводимо 


результати лише для /-компоненти). Для порівняльного аналізу побудовано 
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аналогічні карти за результатами розрахунків, зроблених на підставі 
міжнародної аналітичної моделі ІСВКЕТО 251. 

Побудовані карти ізопор дали змогу дослідити особливості просторової 
структури вікового ходу 1 Його динаміку на зламі ХХ-ХХІ ст. На рис. 1-4 
зображено карти ізопор /-компоненти для послідовних 5-річних часових 
інтервалів. У структурі  5У  7/-компоненти | в 1985-1990 рр. за 
експериментальними даними 1 за даними, розрахованими на підставі моделі 
ІСВКЕЇО, чітко виділяється Балканський фокус ЗУ, інтенсивність якого 
становить 24 нТл/рік, а єпщентр є в районі Аппенін-Балкан (див. рис. 1). Така 
просторова структура 5У практично без змін збереглась 1 для наступної 
епохи 1990-1995 рр. (див. рис. 2). Зазначимо про хорошу узгодженість для 
цих періодів модельного поля 5У з даними магнітних обсерваторій. 

Різкі зміни в структурі 5У сталися в 1995-2000 рр. ( див. рис. 3). У цей 
період відбувся розпад Балканського фокуса 5У. Вже в 2000 - 2005 рр. на 
території Європи спостережено доволі спокійну структуру 5У з монотонним 
збільшенням вікового ходу в північно-східному напрямі (див. рис. 4). Є всі 
підстави припустити, що у північно-східній частині Європи відбувається 


формування нового фокуса 5У. 



































те 





Рис. 1. Віковий хід /-компоненти геомагнітного поля в 1985-1990 рр. Ліворуч 


- дані моделі ІСВЕ-10, праворуч - дані магнітних обсерваторій; значення 


ізопор, нТл/рік. 
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Рис. 2. Віковий хід 7/-компоненти геомагнітного поля в 1990-1995 рр. 


Позначення ті ж, що й на рис. І. 
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Рис. 3. Віковий хід 7-компоненти геомагнітного поля в 1995-2000 рр. 


Позначення ті ж, що й на рис. І. 
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Рис. 4. Віковий хід /-компоненти геомагнітного поля в 2000-2005 рр. 


Позначення ті ж, що й на рис. І. 


Отже, у структурі вікового ходу /-компоненти магнітного поля Землі на 
території Європи виявлено фокус 5У регіонального масштабу, який розпався 
під кінець ХХ ст. Час його існування не перевищував 25-30 років. Відповідно 
до просторових 1 часових характеристик, цей фокус належить до 
короткотривалих 1 пов'язаний з генерацією магнітного поля в земному ядрі. 
Просторова структура 5У /-компоненти на території Європи в загальних 
рисах досить добре апроксимована моделлю ІСЕКЕ-10, водночас для окремих 
компонент, особливо для магнітного схилення 0), є досить суттєві відмінності 
даних магнітних обсерваторій від результатів, розрахованих з використанням 
згаданої моделі. 

Робота виконана за сприяння Державного фонду фундаментальних 
досліджень України в рамках наукової теми "Вікові варіації геомагнітного 
поля в Європі" Мо Ф28. 6/043. 
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РЕСПЇ.АВІТІЕ8 ОЕ 5РАСЕ-ТІМЕ 5ТЕОСТОВЕ ОЕ АСЕ УАВІАТІОК ОК 
ЕАВТН МАСХЕТІС ЕІКІ.Д ОУЕВ ЕОБКОРЕ АТ ТНЕ ЕХР ОК 20" - 
ОМ ТНЕ ВЕСІМХІХС ОК 21!" СЕХТОВІЕ5 
У. Макзітсепик, Хи.Ногодуз Ку, Р. МагсрепКо 


Маря ої аєе 5рійшпє ої єеотаєпейс Пе! омег Бигоре Їог Їоиг 5-уеаг іпіегуаїя 
пот 1985 го 2005 уеаг Брауе Бееп Би. Тре тар5 ої /-сопаропепі аге ргезепіед. ТРе 
їоси5 ої аєе 5пййшпе ої /-сопаропепі, утПпісП уга8 БгоКеп аг Ше епдої 208 сепіиту, Па5 
Ббееп їоппа. П ра8 Ббееп сопсТидеа Шаг Фе 5оигсе ої із Госи5 15 тп Пашіа писівеаг ої Фе 


Кагф. 


МАТРИЧНИЙ МЕТОД У ЗАДАЧАХ ВИЗНАЧЕННЯ МЕХАНІЗМУ 
ВОГНИЩА ЗЕМЛЕТРУСУ 


Д. Малицький 


Карпатське відділення Інституту геофізики ім.С.І.Субботіна НАН України 
Е-паї!: дапупто(Фсб-ієрі.Туту. ма 


Дослідження вогнищ землетрусів є одним з найважливіших завдань, ЯКІ 
виконує сейсмологія. З погляду математичного моделювання динамічні 
задачі механіки руйнувань виявилися нерозв'язні класичними методами. 
Тому постало завдання створити нові підходи, у рамках яких можна 
розглядати явища поширення розривів у неоднорідних середовищах. Досі 
тривають дискусії про точкові моделі вогнища. Розвиток обчислювальних 
методів привів до вдосконалення та узагальнення результатів, що дало змогу 
визначати параметри джерела, використовуючи метод підбору для 
багаторазового розв'язування прямої задачі. Використання функцій Гріна 
допомагає отримати вираз для поля зміщень у будь-якій точці через вектор 


стрибка переміщень |1|. Ця обставина привела до появи праць, у яких 
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аналізували  хвильове поле для окремих випадків розподілу стрибка 
переміщень |2|. Зарубіжні та вітчизняні вчені запропонували низку підходів 
для визначення параметрів джерела |Ї, 2, 4|. Визначення механізму вогнища 
дає відповідь на знаходження нодальних площин, одна з яких збігається з 
площиною розриву у вогнищі. Це відкриває можливості для виконання інших 
сейсмологічних завдань, зокрема побудови моделей підготовки землетрусів. 
Сьогодні механізм вогнища визначають з використанням великої 
кількості станцій для точкового джерела. Однак точкова модель джерела як у 
просторі, так 1 в часі не дає відповіді на важливе питання - наприклад, 
визначення часової функції джерела (5ТЕ). Побудова механізму вогнища 1 
визначення сейсмічного тензора як функції часу є нашою метою. Для 
виконання завдання | використано - матричний метод, щоб виконати 
математичне моделювання хвильових процесів у шаруватих середовищах. 
Розглянемо поширення сейсмічних хвиль у шаруватому середовищі 
зумовлене точковим джерелом у вигляді тензора сейсмічного моменту. Поле 
переміщень на поверхні в дальній зоні визначене з використанням 


матричного методу |-3 |: 


и) 3 о ко сої М 0 б о при М 
КО) -У|к ЦІ. М геї АК, по - Їх ЛІ У іо Ак. (1) 
Мро/ із10 і-50 
де 

ло 90 72, Ш 0 

пе) ко Жені , 13-15; 

оре 5 іс 8 «9, Дал, 
біг 


"9, Л - функції Бесселя; К - горизонтальна компонента хвильового числа; 
Фа - азимут; 


Му М у С05Фа М у 8 Фу, Мо ЗМ іо, 


М з -со8"Фи Мові? Фа М уу 50 2Фи М ду, 
М4 з -с052Фа М п С052Фи М у 25 2Ф М у, 


М5 з М у СО5Фи М 8 Фа, 
Ме з-8іп 2Фа М у 8 2Ф М у о 2С0852Фи М у; 
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коефіцієнти 2 визначені з характеристичної матриці середовища |5|; М; - 
компоненти тензора сейсмічного моменту. 
З використанням аналітично-числового підходу |6| тензор сейсмічного 


моменту визначений за даними записів на станції у вигляді 
Мак, (2) 


де матриця К визначена швидкісною моделлю середовища; 00) - вектор 


переміщень, зумовлений прямими Р 1 5 хвилями на поверхні середовища. 
Визначення для тензора М, згідно з (2), має єдиний розв'язок як розв'язок 
системи з шести рівнянь з шістьма невідомими. На практиці матриця К є 
погано зумовленою. Тому використання сейсмічних даних для події, яку 
зареєструвала більша кількість станцій, приводить до надійнішого результату 
для сейсмічного тензора. 

Нехай для певної події є У станцій, на яких записані сейсмограми для 


конкретного землетрусу. Тоді для кожної станції запишемо матричне 


рівняння: 
- 110) 
КкмМ-Юе, 
- 110) 
Км й оФ 
МО" УбУ 


Система матричних рівнянь (3) є перевизначеною щодо шести 


компонент тензора М, Її перепишемо у вигляді 

смород, (4) 
де матриця С- Ж, -- К, / розмірністю 6х6М, тому що кожна з матриць К; 
(п «1,М) має розмірність 6хб. Вектор у розмірністю І1х6М, оскільки 


вектор бо (п - 1, М) має розмірність Іхб, 


Перевизначена система може бути несумісною, тому зведемо систему 


лінійних рівнянь (4) до простішого вигляду з квадратною невиродженою 


матрицею СС, помноживши (4) зліва та справа на матрицю С. яка щодо 
матриці Є є комплексно-спряженою та транспонованою: 
с'єм-в ою. (5) 
або 
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м-«С сус об (6) 


Матричне рівняння (5) називають нормальним рівнянням | 1 |. Розв'язок 
(6) для тензора М названо розв'язком узагальненої оберненої задачі. Отже 
отримано розв'язок для визначення тензора сейсмічного моменту (6), для 


якого серед усіх можливих розв'язків мінімізується норма 
р 
о -смі т 


На підставі (6) 1 (7) доходимо висновку, що розв'язок узагальненої 
оберненої задачі відповідає методу найменших квадратів. Зазначимо, що 
розроблений метод є ефективним для слабких місцевих землетрусів, хоча 
його використання потрібно перевіряти як на тестових прикладах, так 1 на 
реальних записах. Крім того, розв'язок (6) дає змогу визначати компоненти 
тензора сейсмічного моменту як функції часу, тобто в задачі розглядають 
точкове джерело за просторовими координатами, проте розподіленими в часі. 

Механізм вогнища. Визначення сейсмічного тензора Й механізму 
вогнища за розробленою методикою розглянемо на тестовому прикладі. 
Нехай задано модель середовища: два шари на півпросторі (див. табл.), для 
якої джерело розміщене на глибині /,-2000 м (параметри другого 1 третього 


шарів однакові). Тензор сейсмічного моменту заданий у вигляді 


5.0ж1019  -14ж1019 -70ж1013 
М «| -14ж1018 д452ж1018  -83хж10ї |. (8) 
-7,021019 -83х1019 -952хж1013 


Часова часова функція 5ТЕ(0) як функція часу в джерелі зобр. на рис. І. 





Рис.1. Часова функція джерела 5ТЕ(б) у вигляді загасальної синусоїди 
Відповідно, кожна компонента сейсмічного тензора має вигляд, як 


показано на рис. 2, а. 





а 6 


Рис. 2. Компоненти тензора сейсмічного моменту як функції часу: а - пряма 


задача; б - обернена задача. 
З рис. 2, а визначено час, протягом якого відбувалася подія: І,-0,7 с. 
Нехай є шість станцій, які записали подію. На рис. 3 показано розміщення 
сейсмічних станцій щодо землетрусу на різних епіцентральних відстанях. 


Нехай модель середовища однакова для всіх шести станцій. 


поря щадні два шари на півпросторі 
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Рис. 3. Розташування шести сейсмічних станцій на відстанях 20, 30, 40, 


50, 60, 70 км від епіцентру. 


На рис. 4, а показані синтетичні сейсмограми в декартовій системі 
координат, отримані внаслідок розв'язування прямої задачі (1), 1 механізм 
вогнища для тензора (5). Унаслідок розв'язування оберненої задачі за 
запропонованою методикою з використанням співвідношень (6-7) визначено 
компоненти тензора сейсмічного моменту як функції часу, які зображені на 
рис. 2, б, а також час тривалості події в джерелі. На рис. 4, б побудовані 
сейсмограми на шести станціях 1 механізм вогнища з використанням 
розв'язку оберненої задачі для тензора на рис. 2, 6. Порівнюючи механізм 
вогнища, отриманий за результатами розв'язання оберненої задачі (див. 
рис. 4, б) з вихідним механізмом вогнища (див. рис. 4, а) доходимо висновку, 
що запропонована методика відображає процеси, які відбуваються в джерелі, 
зокрема, визначено час тривалості події. Зазначимо, що використання 
більшої кількості станцій дає надійніший результат. Крім того, важливо мати 


добрі записи та швидкісну модель середовища. 




















Рис. 4. Синтетичні сейсмограми на шести станціях 1 механізм вогнища за 


результатами: а - прямої задачі, б - оберненої задачі. 
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Отже, розроблено методику визначення механізму вогнища за 
результатами розв'язування оберненої задачі з використанням матричного 
методу. Використання більшої кількості станцій дає надійніший результат 
для визначення тензора сейсмічного моменту. Запропоновану методику 
використовують для слабких сейсмічних подій 1 для невеликих 
епіцентральних відстаней. Інші випадки потрібно додатково досліджувати на 
реальних записах. З використанням розв'язку для сейсмічного тензора 


визначено часову функцію джерела 1 час тривалості події у вогнищі. 
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ОК МЕСНАХІ5ЗМ ОК КОСО5 ЕАВТНОАККЕЗ5 


О. Маукзкуу 
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АНАЛІЗ МІНЛИВОСТІ ГЕОФІЗИЧНИХ ПОЛІВ 
У ГЕОЛОГІЧНОМУ РОЗРІЗІ ЦЕНТРАЛЬНОЇ ЧАСТИНИ 
ПЕРЕДГІР'Я УКРАЇНСЬКИХ КАРПАТ 


Р. Морошан, Х. Заяць 


Львівське відділення УкрДГРІ, вул. Пасічна, 38" м. Львів, Україна 
Е-таї!: паогозрап (ФаКта?ті Їуту ца 


Дентральна частина передгір'я Українських Карпат, через яку проходить 
побудований нами (рис. 1) регіональний сейсмічний профіль МСГТ РП-3 за 
напрямом Перегінськ- -Калуш- Галич, перетинає з південного сходу на північ- 
ний захід частково Скибовий, Бориславсько-Покутський та Самбірський 
покриви Українських Карпат, Косівсько- Угерську підзону Більче-Волицької 
зони та Волино-Подільське закінчення Східноєвропейської платформи. 

Скибовий та Бориславсько-Покутський покриви -- це система структур, 
які нагромаджені одна на одну у формі багатоярусної споруди. Самбірський 
покрив насунутий на автохтонні утворення Більче-Волицької зони. 

Розмита поверхня автохтонних мезозойсько-палеозойських відкладів 
монотонно занурюється в південно-західному напрямі. Регіональні порушен- 
ня у фундаменті (Калуський, Краковецький, Передкарпатський) розбивають 
його на блоки, які сходинкоподібно опускаються під тектонічні покриви 
алохтонної частини Карпат. На схилах виступів ерозійної поверхні спосте- 
рігають перспективні антиклінальні та напівантиклінальні структури у 
міоценових відкладах. 

У межах досліджуваної ділянки передгір'я Українських Карпат у минулі 
роки проведені Західноукраїнською геолого розвідувальною експедицією 
Державної геологічної служби ("Укргеофізика") регіональні геофізичні 
дослідження: сейсмометрії КМЗХ-ГСЗ (профілі РП-ХХПІ, П-ПІ, РП-ХХП, 
СГ-1-67) та МСГТ (регіональний профіль Р-3 Перегінськ-Калуш- Галич), 
гравіметрії, магнітометрії та електрометрії МТЗ. 

Матеріали КМ3ЗХ та МТЗ дали змогу у відділі геофізичних досліджень 
ЛВ УкрДГРІ побудувати картосхему ефективних швидкостей (У), та 
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інтервальних питомих електричних опорів (р) до глибини мінус 5 000 м 


(див. рис. 1). 





з | , | М Ж 
32 Я рено Бп-48 


Рис.1. Картосхема ефективних швидкостей (3) та інтервальних 
питомих електричних опорів (4) центральної частини передгір'я 
Українських Карпат 


1 - границі структурно-тектонічних зон: І - Східноєвропейська платформа, 
ІІ - Більче-Волицька зона, ІТІ - Самбірський покрив, ІМ - Бориславсько- 
Покутський покрив, У - Скибовий покрив; 2 - геотраверс РП-3, 


5 - свердловини глибокого буріння 


Ізолінії ефективних швидкостей 3 (3 300-4 200 м/с) та р 4 (15-20 Ом':м) 
змінюються в південно-західному напрямі згідно з простяганням меж 
структурно-тектонічних зон. 

Електричний опір (рис. 2) закономірно зменшується від 150 до 50 Ом:-м 
у межах Східноєвропейської платформи, досягаючи мінімуму 10-40 Ом:-м у 
Передкарпатському прогині, 1 збільшується до 75 Ом-м у межах покривів 
Карпат. Ізолінії однакових значень ефективних швидкостей зменшуються від 


4300 до 4600 м/с в межах Східноєвропейської платформи 1 до 3 600 м/с у 
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Передкарпатському прогині. На переході до покривів Карпат значення 
ефективний швидкостей знову збільшуються до 4 600 м/с. 

Сумарний аналіз картосхеми на зрізі мінус 5000 м та одержаної нами 
глибинної моделі (до 9000 м) геологічного розрізу Р-3 (див. рис. 2) дає змогу 
простежити кореляційні взаємозв'язки геологічної будови та мінливості 


геофізичних полів. 


Борманавсько-Покувськмії Самбірський покрив Волицька східноєвропейська платфори: 





Перегінськ Львівський палеовойський прогин 
- б -- Грабівка Калуш Галич 
«ні І-Й то В ТТ Р Р сн ай о При Примкі Гролфбрюф 0 Гунфо 0 Грей Грнлюо зара пани плов Пк (км) 
зали 
й кан ва ув вн інн 
і авьеааиь 









і п чсвацй 
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| у 
ре 27 / / 9 787 й 
М. ; дже зн яеєв у дн М 
У: КР хана, / 
ФС 27 ДЕ о ' 7-7 Й 
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Рис.2. Геолого-геофізичний розріз геотраверсу Р-З Перегінськ-Калуш-Галич 


(а) Са 1.Сз РУ «ЕВ 212 ве ОО 
1-ізолінії Уеф в м/с за даними КМЗХ, 2-ізолінії о в Омм за даними МТЗ, 3-графік да в мГл,4-відбиваючі сейсмічні горизонти за даними МСГТ, 


5-заломлена границя Кі поверхні кристалічного фундаменту за даними КМ3ЗХ, 6-тектонічні порушення, 7-фронти покривів та лусок, 
8-ерозійна доміоценова поверхня, 9-можливі глибинні інтрузії 


Ж. КУ У 7ФЯ 72 
у 
г 


ям) 


Тенденції розподілу геофізичних параметрів ефективних швидкостей за 
даними КМ3Х (м/с), інтервальних опорів за даними МТЗ (Ом:м) та графіка 
різницевих аномалій поля сили тяжіння (мГл) до глибини 3 км (див. рис. 2) 
узгоджуються з поширенням у геологічному розрізі літологічних комплексів 
та напрямами смуг структурно-тектонічних зон. Спостережений мінімум 
значень ефективних швидкостей та аномалій сили тяжіння відповідають 
межам залягання в геологічному розрізі міоценових відкладів, перспективи 
нафтогазоносності яких доведені бурінням. Мінімум інтервальних опорів 
охоплює, крім міоценового прогину, палеозойські відклади 
Східноєвропейської платформи, наявність яких у платформному заляганні за 
одержаною моделлю кгеотраверсу Р-3 (див. рис.2) ми прогнозуємо в 
геологічному розрізі автохтону Більче-Волицької зони. 
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З  мінімумами потенційних полів пов'язані зони  розущільнення 
геологічного розрізу, які можна зачислити до перспективних на пошуки 
нафтогазоносних об'єктів. 

Залучення інформації про мінливість геофізичних полів у геологічному 
розрізі дає змогу з урахуванням структурного фактора виділяти у Передгір'ї 
Українських Карпат перспективні на пошуки вуглеводнів смуги та 


прогнозувати напрями геолого розвідувальних робіт. 


АХАТ. У 515 ОК СОРНУЗІСАЇ, ЕТЕТ.О5 ІХ СЕОГОСІТСАЇ, 5ЕСТІТОХ ОК СЕХТКАЇ, 
РАКТ ОЕ КООТНІ, ОК ОККАТХІАХ САКРАТНІАХМ5 


М. Могозбрап, КП. /ауаїз!" 


Гаїогтапоп ароці георпувісаї Пе!48 тп єеоіоєлсаї! 5еспоп аПомуз їо Фф5йпеціяП Фе 
пуадгосагроп58 ргозресйуе 7опез аї Тоо ій ої ОКташіап Саграфіапзя апа го їогесазі Ше 


Фтеспоп5 ої єеоіовісаї! ехріогапоп. 


СТАН І ДИНАМІКА МАГНІТНОГО ПОЛЯ ЗЕМЛІ 
НА ТЕРИТОРІЇ УКРАЇНИ 


М. Орлюк, Т. Сумарук, А. Роменець 


Інститут геофізики ім..С.І.Субботіна НАН України, Київ 
Е-таї!: огіуик Фіорі.Кісу ма 


Квазістаціонарне магнітне поле Землі (МПЗ) є одним із визначальних 
геофізичних чинників щодо процесів у різних оболонках та планети загалом. 
Багато явищ у магнітосфері, іоносфері, атмосфері та біосфері зумовлено 
станом магнітного поля та геомагнітною активністю, тому їхня оцінка та 
прогноз є надзвичайно актуальними. 

Просторово-часова структура індукції магнітного поля Землі В 


визначена сумою полів від різних джерел: 


зд 


ВеВв-В.чвВ. дв, 
де В, -- нормальне (головне) поле Землі, що генерується процесами в 
рідкому ядрі та на межі з мантією 1 яке визначає глобальну просторову та 
часову структуру поля планети; В, -- аномальне магнітне поле (поле 
літосфери), зумовлене, головно, намагніченістю порід; В. -- зовнішнє поле, 
зумовлене впливом сонячного та космічного випромінюваня 1 магнітних 
полів Сонця та навколоземного простору; 08, -- поле варіацій. 

Просторова структура аномалій модуля індукції магнітного поля ДВ 
території України, за даними наземних та аерознімань, має надзвичайно 
складний характер. Розміри та інтенсивність аномалій залежать від 
геологічної будови земної кори та наявності магнітних джерел, зумовлених 
зазвичай залізистими мінералами. 

Безперервну реєстрацію індукції МПЗ та її північної Вх, східної Ву 1 
вертикальної 82 компонент виконують геомагнітні обсерваторіями (ГО). На 
території України тривалий час функціонує три ГО: "Київ", "Львів" та 
"Одеса". За даними спостережень та інтерпретації з'ясовано, що найбільш 
суттєво змінилося головне МІПЗ, зокрема, для території України в 
середньому Вн -1200 нТл за останні 50 років. Це поле надалі ми використали 
для оцінки підмагнічувального ефекту в вікові варіації магнітного поля (5У). 

Розглянемо довготривалі зміни поля вікових варіацій (5У). Величину 5М 
визначають як різницю між послідовними середньорічними значеннями 
трьох ортогональних компонент та повної напруженості геомагнітного поля, 
що безперервно спостерігають на магнітних обсерваторіях. У разі розрахунку 
таких різниць допускають, що варіації поля з періодами менше одного року, 
генеровані зовнішніми джерелами 1 є знакозмінними, автоматично 
виключені. Однак варіації геомагнітного поля, генерованого кільцевим 
магнітосферним струмом, які відображені в горизонтальній (ВН) та 
вертикальній (072) складових поля, завжди одного знака. Тому в 
магнітоактивні роки ВН компонента завжди менша, а Вх компонента -- 
завжди більша, ніж у спокійні роки. Такий же ефект дає сонячно-добова Бз- 
варіація. Внаслідок впливу зовнішніх джерел вікові варіації 5У, отримані за 
середньорічними значеннями, будуть уміщати компоненти, пов'язані з 
сонячною 1, відповідно, магнітною активностями. 

Змінне магнітне поле від зовнішніх джерел індукує в провідних шарах 


Землі струми, магнітне поле яких накладається на дипольне поле. Отже, 


5 


середньорічні значення напруженості геомагнітного поля та його компонент, 
окрім дипольного поля, містять складові від зовнішніх джерел, а також від 
величини "підмагнічувального" ефекту та індукційних струмів у провідних 
шарах Землі. Розділити поле на складові від кожного джерела - надзвичайно 
актуальне завдання. 

Геомагнітні обсерваторії України "Київ", "Львів", "Одеса" розміщені в 
регіонах з різною геологічною будовою, що відображене в інтенсивності 
магнітного поля. ГО "Київ" розташована в районі слабко інтенсивної 
Київської регіональної магнітної аномалії ДАВ-95 нТл, ГО "Львів" - на 
периферії інтенсивної Львівської аномалії ДВ-250 нТл, а ГО "Одеса"? - у 
градієнтній зоні Одеської аномалії ДВ--20 нТл. Отже, треба очікувати різних 
значень "підмагнічувального" та індукційного полів у підстильних шарах 
обсерваторій. 

Для виділення компонент, пов'язаних із зовнішніми та внутрішніми 
джерелами, потрібно від середньорічного значення поля на конкретній 
обсерваторії відняти його дипольну частину, генеровану струмами на межі 
ядро-мантія. Значення дипольного поля В/СКЕ, а також усіх компонент МПЗ 
обчислено за моделлю ЇСВКЕ 1945-2010 Моаєі СоеїПсівпізя 2010. 

Одержану різницю потрібно розділити на дві частини: одна - Ввід 
зовнішніх  магнітосферних та о іоносферних струмів, друга - від 
підмагнічування магнітних джерел 1 струмів у підстильних шарах. Вплив 
зовнішніх джерел на середньорічні значення геомагнітного поля найліпше 
відображений у горизонтальній 8/-компоненті, вплив внутрішніх джерел - у 
вертикальній В/ та повній В компонентах поля. 

Внесок зовнішніх джерел для всіх обсерваторій оцінено за допомогою 
різниці середньорічного значення за всі дні та за спокійними днями (ВН- 
ВЗ). 

Сьогодні для ГО "Київ", "Львів" та "Одеса" з використанням нового 
критерію оцінки збуреності МІЇЗ розраховано аномалії внаслідок " підмагнічу 
вального" ефекту в аномаліях 5У завдяки зміні поля, що намагнічує породи. 
Де випливає з того, що у разі зумовленості сумарної намагніченості породи 
індуктивною компонентою І; «-уН, співвідношення ДАД-(В-Вп)/2Вп-ДВ/2Вп 
повинне бути сталим для будь-якого часового інтервалу. Розрахунки 
засвідчили, що ця компонента поля за останні 50 років становить: -0,75 нТл 
для ГО "Одеса", 2,25 нТл для ГО "Київ" та 6,3 нТл для ГО "Львів". 
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Отже, оцінено просторово-часову структуру МІЗ на території України 
за останні 50 років, запропоновано методику оцінки внеску зовнішніх та 
внутрішніх джерел в аномалії вікового ходу та отримано перші числові 


оцінки за даними українських обсерваторій. 


ЗТАТЕ АХОЬ РУХАМІС5 ОЕ МАСХЕТІС ЕТЕІ.р ОК ЕАДВТН 
ОХ ТНЕ ТЕКІТОБКУ ОК ОККАГХЕ 


М. Огіумк, Т. 5Зипагтик, А. Копіапеїя! 


урасе апа ште 58ігисіите ої дца81-5:айопагу пзаєпейс Пе! ої Багії оп Ше Іеггіїогу 
ої ОКташте дшгіпє Їа5і: 50 уеаг5 Ба5 Бееп ез5итагед. Тре теоай ої езітайоп ої 
сопілібицоп ої ехіегпа! апа 1піегпаї! 5оипгсез8 їпіо апотаПез ої аєе 5рІй 15 ргоро8ей, Фе 
питегіс езштайоп5 от ОКташтіапя обяегуаюогіез аге обіатеай Гог Фе Пг5і те. 
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МЕТОДИ МОДЕЛЮВАННЯ ПРОСТОРОВОГО РОЗПОДІЛУ 
ГОЛОВНОГО ГЕОМАГНІТНОГО ПОЛЯ ТА ЙОГО ВАРІАЦІЙ 


О. Павлюк 


Львівський національний університет імені Івана Франка, м. Львів, Україна 
Етаї!: оКк5 рауїуик Ф уапоо.сот 


Для опису спостережуваного геомагнітного поля та його варіацій 
пропонують різні моделі глобального та регіонального масштабів. Серед 
глобальних моделей часто використовують Міжнародну модель головного 
геомагнітного поля (ІСКЕ), яку на кожні 5 років розробляє Міжнародна 
асоціація магнетизму й аерономії, або Світову модель магнітного поля 
(УУММУ, розробкою якої займаються спільно Національний геофізичний 
центр даних (МОРС, Боулдер Со, США) та Британська геологічна служба 
(ВС5, Единбург, Шотландія) під керівництвом силових відомств цих країн. З 


математичного погляду глобальні моделі є  сферично-гармонічним 
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зображенням головного магнітного поля Землі. Сьогодні в моделі ІСВКЕ-ІЇ 
розраховані коефіцієнти розкладу до 13 порядку |11. 

Мета сучасного регіонального моделювання --- опис геомагнітного поля 
на ділянці земної поверхні, який забезпечує ліпше просторове розділення 
локального поля на окремій території з високою густотою даних. 

Серед регіональних методів моделювання у літературі згадують такі |2- 
6): сферично-гармонічний аналіз, поліноміальний, прямокутний гармонічний 
аналіз, сферичний кеп-гармонічний аналіз та його модифікації. 

Класичною методикою опису геомагнітного поля в просторі без джерел 
є сферично-гармонічний аналіз, який зображає просторову структуру 
потенціалу геомагнітного поля у вигляді усіченого розкладу за сферичними 
гармоніками ряду Фур'є, залежними від довготи, та поліномами Лежандра, 
залежними від косинуса широти. Однак сферично-гармонічне наближення 
найкорисніше у зображенні глобального геомагнітного поля та його варіацій, 
коли дані спостережень розподілені більш-менш рівномірно по всій сфері. 
Проте на обмеженій території поліноми Лежандра та тригонометричні 
функції не є вже найзручнішими базовими функціями для побудови 
потенціальної моделі |3-35|. Сферично-гармонічна модель не відтворює 
регіональних особливостей поля. Виникає потреба розробити регіональну 
модель геомагнітного поля на заданій території. Особливо інтенсивно такі 
роботи ведуть у Європі, Північній Америці, Японії та на Далекому Сході, де 
добре розвинута мережа спостережень 1 накопичена велика база даних 
вимірювань поля. 

Найпростішими регіональними моделями поля є різні поліноміальні 
розклади по широті 1 довготі. Наприклад |2 |: 

Е(Ф.Ю) з ад атфчаз чн азф' наці захо, (1) 
де Е(Ф,А.) - - загальне позначення будь-якої компоненти магнітного поля; Ф -- 
широта; 2 -- довгота; а; -- невідомі коефіцієнти, які підбирають за 
допомогою методу найменших квадратів на підставі даних спостережень. 

Недоліком таких моделей є те, що змодельовані за їхньою допомогою 
потенціали поля не є розв'язками рівняння Лапласа, а отже, не описують 
фізичної природи поля. Інший недолік той, що отримані вирази не залежать 
від радіальної змінної 7, а тому не дають змоги зробити екстраполяцію на 


довільну висоту над поверхнею. Незважаючи на ці недоліки, поліноміальна 
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модель досить корисна 1 завдяки простоті часто вживана, наприклад, коли 
потрібно відстежити наявність аномалій. 

Протягом декількох останніх десятиліть запропоновано нові регіональні 
методики опису геомагнітного поля, у яких поле відображено градієнтом 
скалярного потенціалу. ШПрикладами таких методик є прямокутний 
гармонічний аналіз |6|, у якому використовують прямокутну систему 
координат замість сферичної, та сферичний  кеп-гармонічний аналіз 
2-6), головна ідея якого -- це розклад магнітного потенціалу в змінних 
сферичної системи координат у два ряди за базовими функціями, які 
задовольняють рівняння Лапласа в межах вибраної ковпакоподібної області, 
виокремленої колом певного радіуса на сфері. 

Кеп-гармонічний аналіз (ЗСНА) уперше запропоновано 1985 р. для 
моделювання поля в районі Північного полюса, тобто регіон моделювання -- 
це ковпак, обмежений тілесним кутом, половина якого становить бо, ІНШИМИ 
словами, Ф змінюється в межах від 0 до 2л, а 2. -- від 0 до бо. У підсумку 
отримали крайову задачу, граничні умови якої підбирають так, щоб 
можливити розклад потенціалу в рівномірно збіжний ряд за базовими 
функціями в межах запропонованої області. 

Загальний розв'язок цієї задачі має вигляд 


К прах К Ке пу ЧІ і за ри 

С(г,0,Ю а Ке 53 уз - (: со5тА Й 8іп ті б (со50). (2) 
Ке0 тей 

Потенціал Ї/ залежить від радіальної змінної 7, доповнення до 


географічної широти (колатітуди) 0 та від довготи 4; Ку -- середній радіус 


Землі; Я ПР -- коефіцієнти розкладу; ди -- приєднані функції 
К 
Лежандра з нецілим степенем п; та цілим порядком т. Дробове значення 


степеня означає, що я" не є більше поліномами зі скінченною лінійною 
К 


комбінацією степенів. Вони утворюють ортогональний базис не у всьому 
просторі, а лише в межах виділеного підпростору. 

Серед очевидних переваг цього методу те, що потенціал поля є 
розв'язком рівняння Лапласа 1 те, що вдається зменшити фальшиві значення 
на краях області. Серед недоліків виділяють |3, 6|: 1) значну громіздкість 
пошуку коефіцієнтів розкладу, 2) через дробове значення п; втрачено зміст 


поняття про мультиполь, 3) показник степеня п, пов'язаний з довжиною 
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хвилі, яку відображає відповідна гармоніка, 1 для малих розмірів ковпака, 
виявляється, довжина хвилі становить близько 2 700 км, тобто сферичні 
гармоніки для малого ковпака не відображають явищ у земному ядрі. 

Ті ж самі методи використовують для аналізу вікових варіацій поля. До 
моделювання вікових варіацій існує два підходи: 1) або моделювати значення 
кожної компоненти поля в задані епохи окремо, а потім різницю між двома 
послідовними значеннями поля ділити на відповідний часовий інтервал, 
2) або моделювати часові похідні від потенціалу. 

Сьогодні ведуть дискусії про доцільність 1 межі використання того чи 
іншого методу, а також пропонують різні модифікації цих методів. З аналізу 
літературних джерел |5-7| можна зробити висновок, що сферичний кеп- 
гармонічний аналіз ліпше описує аномальне геомагнітне поле 1 може бути 
хорошим методом також для головного поля, однак потребує врахування 
більшої кількості коефіцієнтів розкладу, ніж інші регіональні методи, що 
поряд із математичною громіздкістю робить його менш привабливим для 
дослідників. Існує також думка |6|, що можна просто перераховувати 
коефщієнти у звичайній сферично-гармонічній моделі, враховуючи дані для 


заданої території 1 вносячи поправку до ІСВЕ моделі. 
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ЗРАТТАТ, МОРЕШЛІХС МЕТНОРЗ5 ОЕТНЕ МАТХ СКОМАСХЕТІС ЕТЕІ., 
АМО ІТ5 УАВІАТІОМУ5 
О. РаміушКк 
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Розвиток фундаментальних основ фізики Землі та геофізичних моделей 
сприяє застосуванню цілого класу різноманітних прикладних задач, що 
виникають у разі освоєння підземного простору. Прямі обернені методи, які 
грунтуються на виборі області параметрів, уникають обмеження лінеаризації, 
у них є лише одна вимога, щоб можна було оцінити модель на узгодження з 
даними. Важливий аспект фізичних наук - зробити висновки про фізичні 
параметри з даних. Загалом закони фізики забезпечують спосіб розрахунку 
значень даних для заданої моделі. Це - пряма задача. В оберненій задачі 
намагаються реконструювати модель з набору вимірювань. В ідеальному 
випадку існує точна теорія, що визначає, як треба опрацьовувати дані для 
того, щоб відтворити модель. Моделювання динамічних процесів, що 
відбуваються в неоднорідному геофізичному середовищі після збудження 
зовнішніми та внутрішніми полями різної фізичної природи (поверхневими 1 


заглибленими сейсмічними джерелами коливань, припливними силами, 
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хвильовими процесами в атмосфері, техногенними процесами тощо), 
пов'язане з потребою врахування всіх особливостей складної будови 
досліджуваних реальних геологічних структур. 

Фізичні та  механіко-математичні моделі жхжвильових процесів у 
неоднорідному  геофізичному середовищі широко застосовують для 
розв'язування прикладних задач геофізики, геомеханіки та сейсмології. 
Неоднорідність 1  несуцільність геологічних середовищ є причиною 
концентрацій амплітуди хвильових полів на малих дефектах структури. 
Сильна нелінійність фізичних властивостей гірських порід виявляється під 
час поширення слабких сигналів - сейсмічних 1 електромагнітних. Набір 
невідомих величин - це модельний вектор у просторі параметрів. 
Здебільшого використовують просту скалярну незбіжність для вимірювання 
розбіжності між спостереженнями 1 теоретичними передбаченнями з моделі. 
Функція неузгодження формує багатовимірний "рельєф", крізь який повинен 
пройти алгоритм. Зазвичай неоднозначність задачі 1 наявність помилок у 
даних означає таке: якщо модель узгоджується з даними, то необмежена 
кількість моделей узгоджуватиметься з цими даними. 

У прямому методі збиранням даних обійти неоднозначність задачі 
можна повністю шляхом пошуку всіх класів моделі, що узгоджуються зі 
спостереженнями. Друга стадія оберненості полягає в отриманні інформації 
про параметри моделі шляхом вивчення цього повного класу потенційних 
розв'язків. Наприклад, цього можна досягти, якщо спробувати визначити 
властивості характеристик, розділених всіма задовільними моделями, ЯКІ 
отримані на стадії пошуку. На практиці більшість задач охоплює велику 
кількість невідомих величин (десятки-тисячі параметрів), 1 навіть у разі 
грубої дискретизації кожної параметричної осі це швидко веде до 
астрономічної кількості реальних моделей, інколи до десяти тисяч. У цих 
випадках неможливо відібрати всі моделі, а тому пошук задовільного опису, 
який грунтується на  випадковій  (стохастичній) вибірці, Є завжди 
ефективнішим. Якщо ціна розрахунку функції неузгодження велика 1 
розмірність простору висока, то дуже важко порушувати питання про аналіз 
похибок, неоднозначності й компромісів між невідомими. Часто це зводиться 
тільки до пошуку найліпших даних, які узгоджують модель та оцінки 


похибок, використовуючи лінеаризацію, щоб найліпше працювати на стадії 
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оберненості "аналізу похибок" (або які гіпотези перевіряти) у випадках, коли 
прийнятними є прямі методи. 

Ніяка обернена задача не є повною без опису класу моделей, що не 
заперечує дані; це називають задачею оцінки. Нелінійність ускладнює задачу 
оцінки. Перша причина цього та, що врахування нелінійної похибки є 
важкою задачею; друга причина - для деяких нелінійних задач модельні 
параметри впливають на метод, яким у моделі задають дані. Усі процеси в 
геосферах, у тому числі в земній корі, можливі завдяки енергії від потоку 
сонячної енергії, роботі приливних сил та власному  гравітаційному 
потенціалу. Планета загалом споживає стільки енергії, скільки зможуть 
освоїти її структурні утворення. Структури геосфер регулюють потоки 
енергії, і цю функцію вони виконують шляхом зміни речовинного складу, 
своєї структури 1 підтримки внутрішнього руху. Ці структури забезпечують 
також динамічну рівновагу між припливом, поглинанням 1 релаксацією 
механічної енергії в системах земної кори. Доцільно виділити проблеми, 
вирішення яких пов'язано з розвитком фундаментальних досліджень земної 
кори |3, 4|. Складові частини цієї проблеми такі: 

- розроблення нового розділу геомеханіки та геодинаміки - розвиток 
повільних рухів у структурах земної кори. Ці рухи супроводжуються 
дезінтеграцією середовища 1  характеризовані | стійкими режимами 
деформування. Суть тут полягає в тім, що природні тверді тіла великих 
розмірів у разі дуже повільного деформування, хоча 1 зберігають локально 
пружні властивості, загалом поводяться як в'язкі рідини; 

- виявлення фізичних механізмів генерування напруженого стану 1 
процеси його релаксації. Природа сил, що формують напружений стан на 
різних масштабних рівнях, може бути різною. Це 1 гравітація, 1 теплова 
конвекція, 1 приливні сили, 1 сили інерції, 1 планетарні причини. Неочевидна 
природа рушійних сил у тектонщці літосферних плит, а також сейсмічний 
режим як наслідок деформаційних процесів у земній корі; 

- вивчення процесів формування структур руйнування, блокової 
структури й головних структурних елементів земної кори: зон розтягання 1 
стискання, підняття Й опускання, тектонічних розломів різного рангу, 
горизонтального розшарування тощо. Вплив деформаційних процесів на рух 


флюїдів по проникних шарах 1 розломах. Поверхневе масоперенесення є 


6Ї 


активним рельєфоутворювальним фактором, особливо в гірських районах, де 
розвиваються зсувні явища й утворюються лавини. 

Фундаментальні основи фізики земних процесів. Проблеми, пов'язані 
з розвитком теоретичної науки 1 прикладної геофізики, відображають фізичні 
уявлення про деформаційні властивості твердих гірських порід (тверде тіло, 
пружне тіло, сипуче середовище, в'язко-пружне середовище, пластичне 
середовище, тріщинувата порода, неоднорідне середовище), записані у 
вигляді відповідних рівнянь стану. Наведені характеристики фізичних 
процесів у структурних оболонках Землі ставлять перед геологами, фізиками 
та геофізиками проблему з вирішення та розробки алгоритмів обернених 
задач для механізмів формування властивостей геологічного середовища. 

Традиційні методи експериментальних геофізичних досліджень, що 
грунтуються на пасивному моніторингу природних полів 1 природних 
збурень, які мають випадковий характер, таких як землетруси, сонячні 
спалахи, атмосферні електричні розряди й інші явища, потребують довгих 
рядів спостережень 1, зазвичай, виявляють лише імовірнісні закономірності. 
Дослідження взаємозв'язку між геофізичними полями різного походження 
має фундаментальне значення як для поглиблення знань про природу, так 1 
для вирішення найважливіших проблем, таких як прогноз землетрусів, 
дослідження геодинамічних процесів. 

З огляду на це принципи фізичного розуміння геологічних процесів та 
створення відповідних фізичних моделей реальних геологічних структурних 
утворень та оболонок Землі полягають у такому: 

е Земля - космічне тіло, 1 процеси у ній підпорядковані фізичним 
законам, що можна формалізувати рівняннями математичної фізики; 

е головні чинники тектонічної активності Землі, зокрема гравітаційна 
диференціація речовини, мають внутрішню природу; 

е Земля та її оболонки (кора, мантія 1 ядро) є спільною фізичною 
системою, а тому необхідно загалом розглядати горизонтальні 1 вертикальні 
зв'язки між ними; 

е геодинаміка Землі є наслідком фізико-хімічних та механічних 
процесів, що відбуваються у надрах Землі; 

е нерівномірна геодинамічна активність у просторі 1 з часом приводить 


до ймовірнісної періодичності головних процесів еволюції планети. 
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Особливості поведінки геофізичних, гідрогеодинамічних 1 геохімічних 
процесів не можна описати в рамках лінійних моделей математичної фізики. 
Нелінійність є проявом сукупності багатьох характеристик 1 стану речовини 
середовища. Тому нелінійні сейсмічні 1 геофізичні процеси використовують у 
прикладній 1 фундаментальній  геофізиці для дослідження процесів 
"самовзаємодії", складних "ланцюгів" взаємодій, що охоплюють різні 
геофізичні поля та процеси в земних оболонках. 

Наслідком такого підходу є поняття про геофізичне середовище, якому 
притаманні: 

е ієрархічна неоднорідність, несуцільність, різнорангова диференціація 
від малих частинок зерна мінералу до планетарних неоднорідностей; 

е фізична нелінійність характеристик 1 процесів; 

е активність, тобто здатність постійно випромінювати енергію у формі 
сейсмічної, акустичної Й електромагнітної емісії або тепла, що можна 
вимпряти експериментально; 

е мінливість фізичних властивостей з часом, зумовлена також 
зазначеними активністю 1 фізичною нелінійністю; 

е здатність геофізичних процесів взаємодіяти між собою безпосередньо, 
наприклад, сейсмоелектричний ефект, та в опосередкованій формі через 
низку складно пов'язаних між собою взаємодій. 

Проблеми вибору обернених задач та методів у геофізиці. Обернений 
процес близько пов'язаний з прямим моделюванням. Пряме моделювання 
грунтується на використанні математичного взаємозв'язку, такого як 
хвильове рівняння, для синтезування відгуку моделі Землі на збудження, 
наприклад, імпульс сейсмічної енергії. Такі моделі визначають у зв'язку з 
поняттями набору параметрів для вибраного шару, таких як шарова 
швидкість та густина шару. Отже, звичайно важливо вибрати процедуру 
прямого моделювання, яка може адекватно описати фізичні процеси 
спостереження. У випадку сейсмічного дослідження пряме моделювання 
пов'язане з алгоритмом, що виробляє синтетичну сейсмограму із незначною 
різницею для скінченного розв'язку хвильового рівняння та вимірюваних 
величин. У разі гравітаційних даних процедура прямого моделювання 
передбачає розрахунок гравітаційних полів за відомим приповерхневим 
розподілом густини досліджуваного шару структурних оболонок Землі. Для 


вибору потрібної математичної моделі також важливо знати, скільки 
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модельних параметрів треба використати 1 які параметри є найвагомішими у 
впливі на розв'язок конкретної фізичної задачі. 

Вибір правильної моделі напряму залежить від цієї задачі розвідки. 
Оберненість такого моделювання спонукає реконструювати характеристики 
приповерхні з заданого набору вимірювань 1 робити це шляхом /" підгонки" зі 
спостереженнями з точністю, не вищою від похибки вимірювань. Ось чому 
вибір "доброї" моделі є ключовим моментом. Навіть якщо припустити, що 
вибір моделі адекватний, все ще є багато проблемних моментів. За останні 
два десятиліття обернені задачі геофізики складалися здебільшого з розробки 
методів лінеаризації нелінійних задач. Лінеаризовані методи багаторазово 
поліпшують початкову гіпотезу в разі невідомих параметрів, якщо 
використовувати локальну похідну інформацію. Часто похідні розрахунки є 
дуже дорогими 1, на жаль, нема гарантії того, що ітерації збіжаться до 
сприйнятливої моделі (тобто набір невідомих величин). Навіть якщо 
збіжність досягнуто, підсумкова модель може бути далекою від правдивої, 
погано узгоджуватися з даними, мати сильну залежність від початкового 
припущення за невідомих параметрів. Якщо рівень нелінійності значний, то 
тоді два або більше цілком різних класів моделі можуть існувати 1 
задовольняти умови. Якщо ж лінеаризаційні методи не спрацьовують, то тоді 
необхідно використовувати глобальніші підходи. Огляд методів теорії 
обернених геофізичних задач для лінеаризованих задач можна знайти в |2, 31. 

Геофізики постійно працювали над розв'язуванням обернених задач від 
початку математичних спроб описати фізичні процеси в Землі. Геофізичну 
інтерпретацію роблять на підставі аналізу спостережуваного набору даних 
про геологічні властивості, таких як сейсмограми або записи потенціальних 
полів. Груба модель процесу, з якого отримують експериментальні дані, 
міститься всередині мозку інтерпретатора; потім інтерпретатор використовує 
цю грубу уявну модель, щоб отримати приповерхневі властивості зі спосте- 
режуваних даних. Тобто висновки про приповерхневі властивості зі спосте- 
режуваних даних є ідентичними з розв'язком так званої оберненої задачі. На 
противагу, пряма задача складається з визначення даних, які отримують з 
відомих параметрів будови приповерхні та за припущень, що виконуються ці 
закони фізики. Дотримуватимемось зазначеного вище широкого означення 
оберненої задачі: тоді всі відомі алгоритми, які можна застосовувати в 


обчислювальному процесі, можна розглядати як процедури інвертування 
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геофізичних даних. Наприклад, дані про поширення сейсмічних хвиль дають 
змогу реконструювати приповерхневу геометрію земної кори | І, 4, 51. 

Розрахунок варіації амплітуди зі зсувом (ВАЗ) | | стосується інверсії для 
властивостей породи, отриманої з вимірювань амплітуди поверхневих 
сейсмічних хвиль. Інверсія може стосуватися різних типів геофізичних 
даних. Можна одночасно або послідовно підганяти набір різних геофізичних 
даних, таких як сейсмічні характеристики середовища, потенціальне поле, 
дані свердловинних досліджень, з такою самою моделлю Землі, й у кожному 
випадку ми припускаємо, що виконується закон фізики. У разі сейсмічної 
інверсії, наприклад, законом є хвильове рівняння або одне з його числових 
наближень. Алгоритми, які грунтуються на фізичному законі, дають нам 
змогу інвертувати спостережувані дані для приповерхневих характеристик. 
Ми подаємо дуже широке означення інверсії, якого нема в загальній 
літературі, тим більше у працях з розвідувальної геофізики. Швсля 
наближеного аналізу теоретичного змісту оберненого методу ми обговоримо 
деякі популярніші процедури в розвідувальній геофізиці, що їх загалом 
ідентифікують як обернені методи. Раніше більшість використовуваних 
геофізичних розрахункових методів розглядали як спроби розв'язати 
обернену задачу: маємо геофізичні дані, абстрактну модель процесу 1 
шукаємо алгоритми, що дають змогу інвертувати її. 

Отже, головні проблеми 1 напрями фундаментальних фізичних 
досліджень в обернених задачах геофізики такі: 

вивчення сучасних рухів 1 розвитку деформаційних процесів у земній 
корі, у тому числі повільних рухів, що супроводжуються нагромадженням 
деформацій 1 концентрацією напружень на неоднорідностях; 

вивчення взаємодії із земною корою потоків енергії, що надходять до неї 
з надр унаслідок розсіяння гравітаційної, кінетичної 1 теплової енергії 
планети; 

вивчення механізмів виникнення 1 перетворення блокової структури та 
функціонального призначення ієрархічних структур земної кори; 

вивчення механізмів формування 1 релаксації тектонічних напружень у 
неоднорідному блоковому середовищі на різних масштабних рівнях; 

розробка геомеханічних 1 геодинамічних моделей структурних оболонок 


Землі, спрямованих на побудову кількісної теорії з урахуванням неоднорідної 
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1 блоково-ієрархічної структури та особливостей розподілу густини речовини 
мантії 1 ядра. 

Теоретичні 1 розрахункові аспекти теорії оберненої задачі корисні для 
розвитку технології опрацювання геофізичних даних про відповідне 
середовище, проте ітеративну геофізичну інверсію наразі ще широко не 
використовують у розвідувальній галузі через надто високу вартість 
отримання спостережуваних даних та складність обчислювальних процедур. 
Адже лише тоді тривимірна глибинна томографія стане економічно 
вигідною, коли алгоритми обернених задач дадуть змогу геофізикам 
отримувати дані спостережень не тільки для приповерхневої структурної 
геометрії кори, а також для щораз більшої кількості необхідних детальних 


фізичних, хімічних та геологічних характеристик. 
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ІМАСК ЕХНАХСЕМЕХТ ОК ІІХЕАМЕХТУЯ 
О5ІХС ЕІСЕХАХАТ. У 515 ОЕ МАСХЕТІС СВАРІЄХТ ТЕХ5ОК 


ВШепі ОКОС 
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Сатриб, 41380 Гг.ті/Косаєїї, ТигКеу 
Е-таї!: риїогФКосаєїї едай.їг 
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